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I. SELF-SIMILAR RENORMALIZATION 



» I ' An approach is suggested for analyzing time series by means of resummation techniques of 

f> |_ theoretical physics. A particular form of such an analysis, based on the algebraic self-similar renor- 

malization, is developed and illustrated by several examples from the stock market time series. 
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Many data in different sciences are presented in the form of time series. The problem of analyzing the latter 
consists in understanding the dynamics of motion from one temporal point to another in the past and, hopefully, in 
forecasting the data for at least some near future. The standard approach to time-series analysis is to try to guess 
what stochastic dynamical system is behind the data, that is, to attempt to model the events producing the considered 
data by stochastic differential or finite-difference equations [1,2]. Such an approach has been proved to be reasonable 
when applied to a system close to a stationary state, but sudden changes in the behavior of the dynamical system 
generating the time series cannot be accurately grasped. 

In this communication we advance a novel approach to analyzing time series, based on resummation techniques that 
are used in theoretical physics. As examples for illustration, we opt for time series generated by economic systems, 
including sharp changes in their behavior. Let us note that different physical analogies and techniques are now often 
7-H ' used for describing ecomonic phenomena [3-8]. 

f***- | As a starting point, we present a discrete set of data, available from the system past, in the form of a polynomial ( 
0^ ■ i.e. as a formal power series), which is a direct analog of a perturbative expansion valid as time ( -> 0. The rest consists 
in the application of one or another resummation technique to the asymptotic expansion where the role of a coupling 
constant is played by time. If only a few points from the system history are known, the degree of the asymptotic 
expansion will be low, and the most popular Pade' summation [9] is very difficult, if possible, to apply. However, 
there is an approach called the self-similar approximation theory [10-17] which successfully works even for a small 
number of asymptotic terms. Here we shall use a variant of this approach, the algebraic self-similar renormalization 
[15-17]. We consider below several examples taken from the history of real stock markets using only a few, up to six, 
O ■ points from the history in order to calculate the following values of the share prices. 

First, let us give the general scheme of the method we suggest. Assume that the values of the sought function, f(t), 
•rH , are known for n equidistant successive moments of time, t = k = 0, l...,n— 1, so that 

H| /(0) = oo, /(l) = ai, /(n-1) = <!„_!. (1) 

_ _ i 

Let us be interested in the value of the function /(n) in the following n-th moment of time. To proceed further, it is 
important to find a compact representation of the set of the data (1) in the form of an explicit function. To this end, 
we can always use a formal Taylor series of the sought function, f(t) = J2k=o A-kt k ■ Since only a finite set of values 
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(1) is available, one can reconstruct only a finite set of coefficients, A/. — Ak(ao, a\, . . . , a„_i), k = 0, 1, . . . , n — 1, 
from condition (1). Thus, we obtain an approximate expansion 



n-l 

f(t) ~ A fe (a ,ai,...,a„_i) t k . (2) 

fc=0 

Let us stress that expansion (2) has no direct sense if continued straightforwardly to the region of finite arbitrary 
t. The problem of reconstructing the value of a function in some distant moment of time from the knowledge of 
a finite set of its values in preceding time moments becomes now equivalent to the reconstruction of the function 
for the finite value of its argument, t = n, from the knowledge of its asymptotic expansion as t — > 0. In theoretical 
physics, such a problem is called renormalization or resummation problem [9,18]. Thus, the problem of forecasting the 
future values from the set of historical data given in the form of time series becomes equivalent to the renormalization 
(resummation) of asymptotic series. An analytical tool for the solution of this problem, called the algebraic self-similar 
renormalization, has been recently developed [15-17]. We describe here only its principal points that are important 
within the context of this letter. The polynomial representation (2) gives for the sought function f(t) the following n 
approximations Pi(t), i = 0,1, ...,n — 1, 

Po(t) = A = OQ, pi(t) = po(t) + Att, Pn-l(t) =p n - 2 (t) + An-it"- 1 . (3) 

The algebraic self-similar renormalization starts from applying to the approximations (3) an algebraic transformation, 
defining Pi(t,s) — t s pi(t), i = 0,1, ...,n — 1, with s > 0. This transformation raises the powers of series (3), and 
allows us to take effectively into consideration more points from the system history. We use in what follows the 
strongest form of such a transformation when s — > oo, which results in a nice exponential representation for the 
sought function. The sequence of so transformed approximations Pi(t,s) is considered as a dynamical system in 
discrete time i — 0, 1, n — 1. In order to describe the system evolution with time, we introduce, according to Refs. 
[10-14], a new variable p and define the so-called expansion function t(ip, s) from the equation Po(t, s) = a t s = ip, 

which results in t(ip, s) = (p/ao) 1 ^- Then we construct an approximation cascade [10-14] whose trajectory points are 
given by the expressions yi(ip, s) = Pi(t(ip, s), s). Embedding this cascade into an approximation flow, one can write 
the evolution equation in the form of the functional self -similarity relation yi +p {ip, s) — yi(y p (<p, s), s). At this stage, 
we can check the effectiveness of the algebraic transformation by calculating the local multipliers, 



7Bj(i, s) 



(4) 

f=P (t.s) 



dip 

as s — > oo. For calculations we use the integral form of the self-similarity relation, 

dip 

where the cascade velocity Vi(ip,s) — yi(ip,s) — yi-\{ip,s) and t is the minimal number of steps of the procedure 
needed to reach the fixed point P*(t, s) of the approximation cascade. It is possible to find P*(t, s) explicitly and to 
perform an inverse algebraic transform after which the limit s — > oo is to be taken. The first step of the self-similar 
renormalization is completed. Then the procedure can be repeated as many times as it is necessary to renormalize 
all polynomials. This is the main idea of the self-similar bootstrap [17]. Accomplishing this program, we come to the 
following sequence of the self-similar exponential approximants 

f*{t,r) =A cxp(^tcxp(^t...cxp(^-Tt^ , j = 1,2,..., n-l. (5) 



A VAi " \A 



In order to check whether this sequence converges, we have to analyze the corresponding mapping multipliers. From 
the equation f*(t,l) = ip, we find t(ip) = {A /Ai) \n(p/A n ). Then we construct an approximation cascade, as is 
described above, getting Zj(ip,r) = f*(t(ip),T), and define 



Mj(t,T) = 



dz 3 (ip,T) 
dp 



(6) 

¥>=A*(t,i) 



Selecting from the sequence of fj two successive terms with smallest Mj(t) = Mj(t, 1), as a rule the two last terms 
of the sequence, we can finally determine t from the minimal velocity condition written as the minimal difference 
condition 
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mm|/;(i,r)-/;_ 1 (t,r)|, (7) 

whose solution r = r(i) allows us to write the final self-similar exponential approximation for the sought function 

f*{t)=f*{t,r{t)). (8) 

The approximants /^_ 2 («, 1) an d ./r*-i( n > 1) usually frame the optimal value /^_i(n). One should analyze, following 
the procedure described above, a set of self-similar approximations for different number of points from the history, 
i.e., one has to calculate several possible f'j{t), j = 2, 3, . . . , n— 1, and to choose among them that which corresponds 
to the smallest mapping multiplier calculated at the fixed point. The value of the multiplier at the fixed point can be 
defined as 

M * (t) = M 3 (t,T(t))+M 3 _ 1 (t,r(t)) ^ (9) 

Now we pass to the illustration of the method by examples from the history of various stock markets. We concentrate 
on different stock market crises, when the prices changed sharply during the period of time comparable to the resolution 
of time series. The data, unless stated otherwise, are taken from the books of International Financial Statistics issued 
by the International Monetary Fund. 



II. MECHANISM OF CRASH 



Consider the behavior of the Dow Jones index in the vicinity of the crisis of October 27, 1997. We intend to give 
an illustration of the self-similar renormalization scheme presented above and also to make some general remarks 
concerning the mechanism of crash. We are going to make a forecast for the Dow Jones index for October 25, 1997 
based on different number of points from the system history. 

Three point forecast. The following values are available in the period of time from September 13, 1997 till 
October 11, 1997 taken with the two-week resolution (data are taken from the Yahoo Finance chart): 

a = 7720 {Sept. 13), di = 7920, a 3 = 8000 {Oct. 11). 

The coefficients of the polynomial (2) are A\ = 260 and A 2 = —60. We find the exponential approximants f* (3, 1) = 
8540.8 and / 2 *(3, 1) = 8120.4 and the mapping multipliers Mi (3) = 1 and M 2 (3) = 0.146. From the minimal difference 
condition min T |/ 2 (3,r) — /*(3, r)|, we obtain r = 0.641428. Finally, the self-similar exponential approximation is 
/ 2 *(3) = 8237 and the modulus of the multiplier (9) is |M 2 *(3)| = 0.485. 

Four— point forecast. The following values are available in the period of time from August 30, 1997 till 
October 11, 1997 taken with the two week resolution (data are taken from the Yahoo Finance chart): 

a = 7690 {Aug. 30), ai = 7720, a 2 = 7920, a 3 = 8000 {Oct. 11). 

The coefficients of the polynomial are A\ = —151.667, A 2 = 230, and A3 = —48.333. Repeating the same procedure 
as above, we find the exponential approximants /| (4, 1) = 7689 and f£ (4, 1) = 7646, with the mapping multipliers 
M 2 (4) = -0.013 and M 3 (4) = 0.046. From the minimal difference condition min T |/ 3 (4,r) - / 2 *(4, r)|, we obtain 
r = 0.6026. The self-similar exponential approximant becomes / 3 (4) = 7674 and the corresponding multiplier 
modulus is I M 3 * (4)| = 0.048. 

Five— point forecast. Consider the historical data for the Dow Jones index in the period of time from 
August 16, 1997 to Octobert 11, 1997. The following data are available: 

a = 7780 (Aug. 16 ), on = 7690, a 2 = 7720, a 3 = 7920, a 4 = 8000 {Oct. 11). 

From condition (1), the coefficients of the polynomial (2), with i = 0, 1, ...,4, are 

A 1 = -48.333, A 2 = -120.833, A 3 = 93.333, A A = -14.167. 

The following three higher-order self-similar exponential approximants, f 2 {t, r), / 3 (£,t), and ft{t,r) can be written, 
so that at t = 5 we have f 2 (5, 1) = 0, / 3 (5, 1) = 7472, and /| (5, 1) = 6113. The corresponding mapping multipliers, 
M 2 (5) = 0, M 3 (5) = 0.319, and M 4 (5) = 13.569, signal that the subsequence /|, / 3 *, /| is unstable. But from the 
minimal difference condition min T 1/4(5, r) — / 3 (5,r)| we can still locate the fixed point and find r = 0.6232. The 
self-similar exponential approximation is /| (5) = 7069, with the multiplier modulus |M|(5)| = 0.163. 
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Six— point forecast. Consider the data for the Dow Jones index in the period of time from August 2, 1997 to 
October 11, 1997: 

a Q = 8170 (Aug. 2), a x = 7780, a 2 = 7690, 

a 3 = 7720, a 4 = 7920, a 5 = 8000 (Oct. 11). 

The coefficients of the polynomial arc A 1 = -771.5, A 2 = 582.917, A 3 = -253.75, A 4 = 57.083, and A 5 = -4.75. 
Following the standard prescriptions, wc can find the exponential approximants /|(6, 1) = 7720 and fg (6, 1) = 7138 
and the mapping multipliers M^(Q) = —0.308 and M$(6) = 0.198. From the minimal difference condition r = 
exp (6A5T /A4) we obtain r = 0.7037. For the self-similar exponential approximation we get f£ (6) = 7328, with the 
multiplier modulus |M*(6)| = 0.087. 

All forecasts, except the three-point one, predict a decrease of the Dow Jones index. Normally, in the absence of a 
dramatic event which can alter drastically the market activity, one should expect that the most stable trajectory is 
to be realized. Expecting the sequence of the multipliers {M*}, j = 2, 3, .., 5, we notice that the four-point forecast 
/g (4) = 7674 is optimal from the viewpoint of stability, that is, it corresponds to the minimal mapping multiplier, and, 
thus, this four-point prediction is to be chosen as the final estimate. The Dow Jones index on October 25, 1997 was 
7715. The percentage error of our forecast equals —0.531%. An upward development, predicted by the three-point 
forecast, appears to be less stable than the decaying trajectories. The level achieved by the market after the bubble 
burst corresponds rather to a correction than to a crash. 

When the normal evolution is disrupted by some unexpected negative dramatic events, such as the crisis in Hong 
Kong, which disturbs the self-similar dynamics, then it may happen that the market temporarily selects not the most 
stable trajectory. In its search for a solution, the market may leave for a while the most stable trajectory and jump 
to a less stable one. Then, it can bounce back to the most stable trajectory or it can fall for a while on the least 
stable trajectory corresponding, in this case, to the five-point forecast. During a short period of time, less than the 
time-series resolution, one can observe a rapid succession of all possible trajectories. After the source of the self- 
similarity violation is lifted, one may expect that the market will return to the most stable trajectory. The motion to 
the least stable trajectory is the most dramatic moment, as is seen from the analysis of the five-point forecast. Since 
both starting polynomial coefficients are negative and \A 2 \ » \A\\, an extremely strong tendency to decay, a sort of 
panic, appears resulting in /|(5, 1) = 0; however the situation is remedied due to the higher-order coefficient A 3 > 0, 
which reflects still existing moderately optimistic views; though such an optimism remains limited by pessimistic views 
represented by the negative A4. The finite value given by f%(5) = 7069 is quite low because of such a combination 
of tendencies. We would like to stress that in the case of the five-point forecast the fixed point determined by the 
locally unstable sequence /|, /| lies very close to the boundary separating stable from unstable trajectories. This is 
why the short-term market changes can be vary rapid, first, to the level determined by /|(6) = 7328, and then to the 
level /g (4) = 7674. We conclude that the following reasons caused the October 1997 fluctuations: (1) the burst of the 
upward bubble accompanied by the appearance of several stable decaying trajectories; (2) the location of one of these 
trajectories, corresponding to strong decay, at the boundary between stable and unstable trajectories; (3) the Hong 
Kong crash that pushed the market from the most stable trajectory with moderate corrections right to this unstable 
trajectory, although shortly after this crisis the market promptly bounced back to the nearby stable trajectories. 



III. MARKET CRISES 



We shall consider several examples of crises happened in different stock markets. In all the cases we have performed 
the same steps as in Sections I and II, but only the optimal forecasts are presented below. It turns out that different 
crises require for their description different number of historical points. Let us start from the crises which can be 
described knowing only three points from the system history: 

(i) Consider the dynamics of the average index of the French industrial share prices from 1986 to 1988, with one 
year time resolution, and make a forecast for 1989. The following historical data are available (the price in 1985 is 
taken for 100, 1985 = 100) : 

a = 153.3 (1986), a x = 177.6 (1987), a 2 = 162.1 (1988). 

From condition (1), the coefficients of polynomial (2), with i = 0, 1, 2, can be found giving A\ = 44.2 and A 2 = —19.9. 
The following exponential approximants (5) can be readily obtained: 

f*(t,r) = a exp ( — rt J , f 2 (t,r) = a exp ( — iexp ( ^rt 
\ao J \ao \Ai 
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from where, at t = 3 and r = 1, we have /j"(3, 1) = 364.078 and /|(3, 1) = 191.8. The mapping multipliers can be 
calculated using formula (6), which yields Mi(3) = 1 and M 2 (3) = —0.048, so that the sequence /*, / 2 is locally 
stable [14]. From the minimal difference condition min r |/ 2 (3,r) — / 1 *(3,r)|, which reads as 

{3A 2 \ 
T = cxp j , 

we find t — 0.505. Finally, the self-similar exponential approximation for the sought share price index is / 2 (3) = 
237.34. The actual value of the index in 1989 was 234.9. The percentage error of our forecast is 1.039%. 

(ii) In 1990 — 1992, the average index of the Denmark shipping share prices had the following values (1990 = 100) 

a = 100 (1990), ai = 100 (1991), a 2 = 92 [1992). 

Let us look for the price in 1993. This case can be considered by analogy with the previous one. The polynomial 
coefficients can be readily found being A\ = 4 and A 2 — —4. For the considered index we have /j"(3, 1) = 112.75 and 
/I (3, 1) = 100.60. The mapping multiplier M 2 (3) = —0.089, hence the sequence of exponential approximants locally 
converges. From the minimal difference condition, we find r = 0.35. Thus, the self-similar exponential approximation 
is /|(3) = 104.289. The actual index in 1993 was 105, so the error is -0.677%. 

(Hi) Consider the average index of the Swiss industrial share prices from 1993 to 1995 (1990 — 100): 

a = 137.4 (93), ai = 159.2 (94), a 2 = 166 (95). 

The corresponding polynomial coefficients are A\ — 29.3 and A 2 = —7.5. Let us make a forecast for 1996. Repeating 
literally the same steps as above, we can easily calculate the exponential approximants /* (3, 1) = 260.508, / 2 (3, 1) = 
184.883, and the mapping multiplier M 2 (3) = 0.076. Then from the minimal difference condition we find r = 0.621. 
So, the self-similar exponential approximation is f 2 (3) = 204.394. The price in 1996 was 206.3, which shows that the 
percentage error of our forecast equals —0.924%. 

We continue with the crises which can be described using four— points from the history: 

(i) Consider the dynamics of the average index of the French industrial share prices in 1970 — 1973, with one year 
resolution (1970 = 100): 

a = 100 ( 70), 01 = 95.8 ( 71), a 2 = 107.4 (72), o 3 = 129.7 ( 73). 

Let us make a forecast for 1974. The polynomial coefficients found from condition (2) are A\ = —13.8, A 2 = 10.45 
and A3 = —0.85. The exponential approximants (5) are 

fi(t,r) = a cxp \ —Tt) , 



f 2 (t, r) = a cxp ( — t exp ( — ^rt 
\a a \A 1 




f*(t,r) = a cxp ^^-texp (^-texp (j^rt 

Therefore, at t = 4, we have A* (4,1) = 57.58, / 2 *(4,1) = 97.366, and / 3 *(4,1) = 93.996. The correspond- 
ing multipliers are Mi(4) = 1, M 2 (4) = —0.166, and M%(A) = —0.087. From the minimal difference condition 
min r |/ 3 * (4, t) - / 2 (4,r)|, which reads 

(AA Z 
t = exp — — T 

V M 

we find r = 0.7767. Thus, the self-similar exponential approximation is /g* (4) = 94.885. Our estimate agrees well with 
the actual price equal to 96.6 in 1973. The percentage error is —1.775%. 

(ii) In 1985 — 1988, the average index of the Netherlands general share prices had the following values: 

a a = 100 (85), ai = 128.7 (86), a 2 = 129.2 (87), a 3 = 119.7 (88). 
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We are going to forecast the price for 1989. The polynomial coefficients are A\ — 48.867, A2 = —23.2, and A3 — 3.033. 
The local multiplier m\(t, s) = 1 + A\ (1 + s) t/Aos, as s — > 00 and t = 4, takes the value mi(4, 00) = 1 + 4Ai/Aq = 
2.995, hence the quality of the sequence of Pi, i = 0,1,..., 3, is not good. In this case, we can consider another 
sequence of approximations, which does not include the constant term. The expansion function now is determined 
from the equation Ait 1+S = ip giving t(tp,s) = (ip/Ai) 1 ^^ 8 ^ . The corresponding local multiplier, defined by the 
formula analogous to (4), that is m 3 .i(t, s) = 1 + A 2 (2 + s) t/A\ (1 + s), as s — > 00, t = 4, equals —0.899 < 1. Thence 
the following self-similar exponential approximants can be obtained (see Ref. [17] for more details): 

f 3 (t,r) = a + A 1 tcxp ( 4^ Tt ) > 



fl(t,r) = a + Aiiexp (^-tcxp (^ T *j 

This leads to the estimates /|(4, 1) = 129.264 and /|(4, 1) = 163.417. The mapping multiplier M 3 (4), calculated from 
the formula analogous to (6), equals 0.613, which tells that the sequence f 3 , fl is stable. From the minimal difference 
condition min T |/|(4,r) — /3*(4,t)| wc find r = 0.6952, which results in the self-similar exponential approximation 
fl (4) = 152.206. Our forecast agrees well with the index in 1989, equal to 151.4, thus the error being 0.532%. 
(Hi) Consider the data for the average index of the Mexican share prices in 1991 — 1994 (1990 = 100) : 

a = 190.1 (91), ai = 291.3 (92), a 2 = 325.6 (93), a 3 = 442.1 (94). 

The polynomial coefficients arc A\ = 184.35, A 2 = —108, and A3 = 24.85. Let us make a forecast for 1995. Since 
mi (4, 00) = 4.879 and 7713^(4, 00) = —1.343, wc will proceed in a close analogy to the case (ii). Repeating the same 
steps as above, we can readily calculate the exponential approximants /| (4, 1) = 260.893 and /| (4, 1) = 480.018 and 
the mapping multiplier M$(4) = 0.13. Then from the minimal difference condition, we get r = 0.584. The self-similar 
exponential approximation is /|(4) = 377.695. The index in 1995 was 389.3; so the percentage error of our forecast 
equals -2.981%. 

Five— point crises: (i) Consider the historical data for the average index of the German industrial share prices 
in 1957— 1961, with one year resolution, and let us make a forecast for 1962. The following data are available 
(1958 = 100) : 

a a = 78 (57), a! = 100 (58), a 2 = 171 (59), a 3 = 272 (1960), a 4 = 282 (61). 

From condition (1), the coefficients of polynomial (2), with i = 0, 1, 4, are A 1 = 16.667, A 2 = -12.75, A 3 = 22.333, 
and A4 = —4.25. The higher-order self-similar exponential approximants are 

f*(t,r) = a exp (^-texp ^^iexp (^rt 




fl(t,r) = a exp ^— texp y— texp y— iexp [^ Tt 

At t = 5, we get / 3 *(5, 1) = 226.885 and /|(5,1) = 199.314. The corresponding mapping multipliers, M 3 (5) = 

1.003 and M 4 (5) = 0.687, signal that the subsequence /|, /| is stable. From the minimal difference condition 
min r |/|(5, t) - / 3 *(5,t)| we have 




The latter equation gives r = 0.577. The self-similar exponential approximation /| (5) = 221.268 agrees well with the 
actual value of the index in 1962 equal to 221. The percentage error is 0.121%. 

(ii) The average index of the Swedish share prices in the period of time from May 1987 to September 1987 had the 
following values (1980 = 100): 

a = 726 , ai = 740 , a 2 = 801 , a 3 = 825, a 4 = 877 . 

What was the value of the index in October 1987 ? This case can be considered by analogy with (i). The polynomial 
coefficients are A\ = -74.75, A 2 = 133.792, A 3 = -51.25, and A 4 = 6.208. The estimates for the index are 
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/*(5, 1) = 632.564 and / 4 *(5, 1) = 710.11. The multipliers are M 3 (5) = 0.861 and M 4 (5) = -0.06, hence we conclude 
that the sequence of the exponential approximants locally converges. From the minimal difference condition, we get 
t = 0.6675. The self-similar exponential approximation /|(5) = 695.732. The actual value of the index in October 
1987 was 697. So that the error of our prediction is —0.182%. 

(Hi) Consider the historical data for the average index of the Israel industrial share prices in 1989 — 1993 (1990 = 
100): 

a = 87 (89), 01 = 100 (90), a 2 = 155 (91), a 3 = 297 (92), a 4 = 418 (93). 

The polynomial coefficients are A\ = 45.25, A 2 = —71.625, A 3 = 45.75, and A4 = —6.375. Let us make a forecast for 
the year 1994. Repeating the same steps as above, we find the approximants f 3 (5, 1) = 569.963 and /| (5, 1) = 146.153 
and the mapping multipliers M 3 (5) = 0.608 and M 4 (5) = 0.004. Then from the minimal difference condition, we get 
t = 0.64. The self-similar exponential approximation is /|(5) = 221.18. The actual value of the index in 1994 was 
257. The percentage error of our forecast equals —13.938%. 

Six— point crises: (i) Consider the historical data for the Standard Statistics index of the New York Stock Exchange 
prices in the period from April 1929 till September 1929, with one month resolution, and let us make a forecast for 
October 1929. The following historical data are available (taken from the League of Nations Statistical Yearbook: 
1926 = 100) : 

a = 193 (Apr. 1929), a x = 193, a 2 = 191, a 3 = 203, a 4 = 210, a 5 = 216 (Sept. 1929). 
From condition (1), the polinomial coefficients, for i = 0, 1, 5, are 

A x = 26.683, A 2 = -49.208, A 3 = 28.333, A 4 = -6.292, A 5 = 0.483. 
Two higher-order self-similar exponential approximants write 

t*u ^ ( Al + f A2 + ( A * + ( Ai * 

f 4 (t,r) — a cxp — texp I — iexp I — — iexp I —r-rt 



«o \Ai \A 2 \A 3 



f 5 (t,r) = ao exp (^texp (^exp ^texp ^texp (J-rt 

The estimates for the index at t = 6 are /J (6,1) = 194.884 and /I (6, 1) = 206.722. The corresponding mapping 
multipliers, M 4 (6) = —0.025 and M 5 (6) = 0.021, signal that the subsequence / 4 , ft, is stable. From the minimal 
difference condition min T |/|(6, r) — / 4 (6, r) one has 




This gives r = 0.7182. The self-similar exponential approximation is (6) = 201.692. The actual value of the index 
in October 1929 was 194. The percentage error equals 3.965%. It is worth noting that our lower bound estimate, 
/! (6, 1), practically coincides with the real value. 

(ii) Consider the behavior of the industrial share prices in the New York Stock Exchange characterized by the 
average Standard&Poor index, in the period of time from the second quarter of 1 986 to the third quarter of 1 987 
(1980=100): 

a = 199.4 (II, 86), on = 198.2, a 2 = 201.3, a 3 = 235.1, a 4 = 253, a 5 = 277.5 (III, 87). 

And let us make a forecast for the fourth quarter of 1987. This case can be considered by analogy with (i). The 
coefficients of the polynomial are 

A l = 52.12, A 2 = -103.717, A 3 = 64.096, A t = -14.883, A b = 1.184. 

For the index we obtain the values /|(6, 1) = 202.109 and / 5 *(6, 1) = 226.696. The mapping multipliers are M 4 (6) = 
—0.01 and M$(6) = 0.012, so we conclude that the sequence of approximants locally converges. From the minimal 
difference condition, we find r = 0.7119. The self-similar exponential approximation becomes f£(6) — 215.293. The 
actual value of the index in the fourth quarter of 1987 was 218.3. The forecast error is 1.377%. 
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(Hi) Consider the behavior of the Standard&Poor index in the period of time from February 1990 to July 
1990 (1985 = 100), 

a = 183.4 (Feb. 90), oi = 188.5, a 2 = 189.2, a 3 = 196.4, a 4 = 202.8, a 5 = 204.9 (Jul. 90). 

Let us make a forecast for August 1990, the time of the Persian Gulf crisis. The coefficients of the polynomial can be 
restored giving 

A x = 19.883, A 2 = -25.158, A 3 = 12.783, A 4 = -2.592, A 5 = 0.183. 

Repeating the same procedure as above, we find the exponential approximants /|(6, 1) = 188.986 and 75(6,1) = 
201.789 and the mapping multipliers M 4 (6) = —0.066 and ^5(6) = 0.043. From the minimal difference condition, we 
obtain r = 0.7327. The self-similar exponential approximation is /s(6) = 197.476. The actual value of the index in 
August 1990 was 188.1. The percentage error of our forecast equals 4.985%, and the lower bound given by / 4 (6, 1) 
practically coincides with the actual value. 

IV. ATTEMPT OF PREDICTION 

Finally, consider the behavior of the Standard&Poor, Dow Jones, and NYSE Composite indices in the period of 
time from April 30, 1997 to September 30, 1997, and let us make a forecast for October 31, 1997. 

For the NYSE Composite index we have (from the NYSE Historical Statistics Archive in Internet): 

a = 416.94 (Apr. 30), ai = 441.78, a 2 = 462.44, 

a 3 = 494.5, a 4 = 470.48, a 5 = 497.23 (Sept. 30). 

The coefficients of the polynomial are 

A x = 104.366, A 2 = -155.195, A 3 = 98.434, A 4 = -24.91, A 5 = 2.145. 

Repeating the same procedure as above, we find the exponential approximants /|(6, 1) = 430.124 and /5*(6, 1) = 
520.109 and the mapping multipliers M 4 (6) = —0.022 and M$(6) — 0.03. From the minimal difference condition we 
obtain r = 0.6974. The self-similar exponential approximation is fi- (6) = 478.855. 

For the Standard&Poor index (data are taken from the DBC Online chart) we have 

a = 783 (Apr. 30), 01 = 845, a 2 = 888, a 3 = 948, a 4 = 900, a 5 = 950 (Sept. 30). 

The coefficients of the polynomial are 

Ai = 222.15, A 2 = -306.292, A 3 = 189.75, A 4 = -47.708, A 5 = 4.1. 

We find the higher exponential approximants / 4 (6, 1) = 818.966 and /j*(6, 1) = 1030 and for the mapping multipliers 
we obtain M 4 (6) = —0.021 and M 5 (6) = 0.028. From the minimal difference condition, we find r = 0.6978. Thus, 
the self-similar exponential approximation is /j* (6) = 935.082. 

Similarly, for the Dow Jones index (data are taken from the DBC Online chart) one has 

a = 6950 (Apr. 30), 01 = 7300, a 2 = 7700, a 3 = 8200, a 4 = 7670, a 5 = 8000 (Sept. 30), 

which gives the polynomial coefficients 

Ai = 1477, A 2 = -2291, A 3 = 1528, A 4 = -399.167, A 5 = 35. 

Following the same steps as above, we find the exponential approximants / 4 (6, 1) = 7108 and /j* (6, 1) = 8400 and 
the mapping multipliers M 4 (6) = —0.025 and M$(Q) = 0.04. From the minimal difference condition we obtain 
t = 0.694088. The self-similar exponential approximation is /g(6) = 7788. 

In conclusion, we have shown that resummation techniques of theoretical physics can be successfully employed for 
analyzing time series. The most convenient such a technique, to our mind, is the algebraic self-similar renormalization 
[15-17]. We have applied this approach to analyzing time series corresponding to different stock markets. From 
our point of view, market crises are somewhat similar to critical phenomena in statistical physics. There exists a 
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temporal region around a crisis, where the behaviour of a market, being intimately connected with its behaviour at 
times relatively distant from visible anomalies [19], begins exhibiting specific features, like log-periodic oscillations 
[3-5]. Such a precrisis region around a phase transition, and precrisis effects are similar to precursor phenonema 
caused by heterophase fluctuations [20]. The point of view that market crises have their origin in the collective 
behaviour of many interacting agents and that the stock market crises are analogous to critical phenomena has also 
been promoted by other researchers [5] . The similarity between the crises of stock markets and the critical phenomena 
of statistical systems makes it possible to apply for their description the same resummation methods, such as the 
algebraic self-similar renormalization [15-17]. Farther development of this approach will be presented in the following 
papers. 
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